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ABSTRACT 



Context. During a solar magnetic field reversal the magnetic dipole moment does not vanish, but migrates between 

poles, in contradiction to the predictions of mean-field dynamo theory. 

Aims. We try to explain this as a consequence of magnetic fluctuations. 

Methods. We exploit the statistics of fluctuations to estimate observable signatures. 

Results. Simple statistical estimates, taken with results from mean-field dynamo theory, suggest that a non-zero dipole 
moment may persist through a global field reversal. 

Conclusions. Fluctuations in the solar magnetic field may play a key role in explaining reversals of the dipolar component 
of the field. 
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1. Introduction 

Magnetic activity is a quasi-regular cyclic process during 
the course of which the solar magnetic field changes its po- 
larity in each nominal 11-year cycle. This cyclic activity is 
usually explained as a manifestation of solar dynamo ac- 
tion somewhere within the solar interior. The Sun and its 
activity are more or less axially symmetric. Of course, the 
distribution of sunspots is not completely symmetric be- 
cause of, e.g., the discrete nature of this tracer of solar ac- 
tivity. However, if distributions of this and other tracers of 
solar activity are averaged over a reasonable time interval, 
they becomes almost axisymmetric, and at any time surface 
deviations from axisymmetry in the form of preferred lon- 
gitudes are about 10% (e.g. Berdyugina et al. 2006). This 
is why a description of the solar activity in terms of ax- 
isymmetric dynamo models seems to be a reasonable step, 
and nonaxisymmetric features can be considered as small 
perturbations. Of course, a deeper understanding of solar 
activity in terms of direct numerical simulations for detailed 
solar magnetic configurations, which automatically include 
small-scale nonaxisymmtric features is a very desirable sub- 
sequent step (see, e.g. Brown et al. 2010, whose model has 
a dipole moment whose evolution does hint at some of the 
desired features). 

The scheme described above looks plausible. However 
it implies that the mean magnetic dipole moment of the 
magnetic field vanishes during the course of each rever- 
sal. The point is that axisymmetric spherical mean-field 
dynamos have two possible directions of their dipole mag- 
netic moment, parallel or anti-parallel to the rotation axis. 
Mean-field solar dynamo models exploit this idea massively. 
Almost all mean-field models assume that the magnetic 
dipole moment vanishes at the instant of field reversal. 
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However, on the other hand observers insist that in re- 
ality the magnitude of the solar magnetic dipole moment 
is reduced at the times of its reversal, but does not van- 
ish exactly. Instead its direction moves continuously from 
one pole to the other on a quite complicated trajectory, 
which varies from one reversal to the other (e.g. Livshits 
& Obridko 2006). The topic was investigated in detail by 
DeRosa et al. (2012) who provided convincing evidence con- 
cerning the migration of the solar magnetic dipole moment 
from one pole to the other during the course of a reversal. 

Of course, it can be claimed that the solar dynamo is 
not fully axisymmetric and that weak deviations from axial 
symmetry in solar hydrodynamics and/or dynamo gener- 
ated solar magnetic field are sufficient to explain the rela- 
tively weak dipole magnetic moment at the times of rever- 
sals. The point however is that the mean-field solar dynamo 
is known to be very robust and it is extremely difficult to 
excite a quite substantial nonaxisymmetric magnetic field 
in a solar type mean-field dynamo. 

A radically different viewpoint would be that that there 
is no need to describe solar magnetic field evolution in the 
framework of mean-field dynamos and that if such models 
have problems, it is time to move to direct numerical sim- 
ulations of the non-averaged induction equation, e.g. along 
the lines of Brown et al. (2010). We feel however that it 
would be much better to attempt to resolve the problem 
rather to ignore it. Mean- field dynamo theory is too use- 
ful a tool in understanding these phenomena to be lightly 
discarded. 

We note that the reversals with nonvanishing mag- 
netic moment touches on one more closely related problem. 
The point is that an inclined magnetic moment creates a 
quadrupole magnetic field mode in addition to the nonax- 
isymmetric one. The appearance of a quadrupole compo- 
nent was convincingly addressed in the context of interest 
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by DeRosa et al. (2012). We note however that dynamo ex- 
citation of a quadrupole dynamo configuration by a spher- 
ical dynamo is much less problematic (see e.g. Moss et al. 
2008) than excitation of nonaxisymmetric configurations. 

The aim of this paper is to give an interpretation of ob- 
servations discussed in DeRosa et al. (2012) in the context 
of nonaxisymmetric dynamo configurations and spherical 
mean-field dynamos. We start by revisiting the problem of 
the existence of nonaxisymmetric solutions in the frame- 
work of mean-field dynamos, and confirm that it struggles 
to explain the phenomenology under discussion. Then we 
suggest a minimal extension of the mean-field theory which 
can explain the phenomenon, by adding explicitly a ran- 
dom (fluctuating) magnetic field component (which is any- 
way implicitly present in any mean-field dynamo) to the 
dynamo generated mean field. We demonstrate that such 
model can explain the nonaxisymmetric magnetic dipole 
directions during the reversal time, within the admissible 
parameter range. DeRosa et al. (2012) provide an accessible 
background for the ideas developed in this paper. 

2. Searching for nonaxisymmetries in solar 
dynamos 

The topic of symmetries of spherical dynamo was exten- 
sively addressed in the early years of dynamo studies and 
the basic message from that epoch still remains valid. In 
particular, Radler (1986) and Moss et al. (1991a) investi- 
gated spherical dynamo models which permit nonaxisym- 
metric solutions, and found that solar- type dynamos pref- 
erentially excite axisymmetric rather than nonaxisymmet- 
ric magnetic fields, and that the axisymmetric solutions are 
very stable with respect to nonaxisymmetric perturbations. 
In fact, we did rerun some similar models and confirmed 
the result. A point to be mentioned in the present con- 
text is that the lifetime of a nonaxisymmetric perturbation 
introduced in a dynamo generated axisymmetric magnetic 
configuration can exceed the reversal time, i.e. if the de- 
sired perturbation arises it can survive through the period 
of reversal. 

Of course, a dynamo can in some situations excite non- 
axisymmetric configurations and efforts were undertaken 
to include such possibilities in the framework of solar dy- 
namos (e.g. Ruzmaikin et al. 1988; Moss 1999). Such at- 
tempts are rather isolated exploratory studies and at the 
moment they do not provide strong motivation for further 
reconsideration of significant nonaxisymmetric field genera- 
tion in more-or-less conventional mean-field solar dynamos. 
We can note that, for example in Moss (1999) the axisym- 
metric and nonaxisymmetric field components oscillate in 
phase. 

Note that the issue of possible excitation of a 
quadrupole axisymmetric magnetic configuration is quite 
different from that of the existence of nonaxisymmetric 
solutions. Of course, the solar dynamo is observed nor- 
mally to excite a dipolar configuration. However a mod- 
est variation in the profiles of the dynamo drivers and/or 
control parameters can result in excitation of a configu- 
ration with quadrupolar symmetry; this problem was re- 
cently revisited by Moss et al. (2008), among many others. 
Moreover, archival solar activity data for the XVIIIth cen- 
tury (Arlt 2009) give a hint that the solar magnetic field 
had quadrupolar symmetry at that time (Illarionov et al. 



2011). A substantial deviation from dipole symmetry with a 
significant admixture of a quadrupole field is known to have 
been present at the end of the Maunder Minimum (Sokoloff 
& Nesme-Ribes 1994). Nevertheless, the above remarks do 
not provide evidence for quadrupole-like type features play- 
ing a role in solar magnetic field reversals. 

We do not give here a detailed verification of the above 
results, but restrict ourselves to just one example. 

For orientation, we performed a simple minded numer- 
ical experiment. We took as a basis the nonlinear nonax- 
isymmetric dynamo code described in the investigation of 
solar active longitudes by Berdyugina et al. (2006). This 
code uses an approximation to the solar rotation law in 
the "convection zone" proper (fractional radius r > 0.7), 
with an isotropic alpha-effect. The model also includes a 
rather thick overshoot layer with slightly reduced diffusiv- 
ity, which allows the angular velocity to become uniform at 
its base (r = 0.64), in a smooth manner. 

First the code was run with slightly supercritical pa- 
rameters, to obtain a stable, oscillating solution with pure 
dipole-like parity. In our dimensionless units, the period of 
oscillation in energy is Pe ~ 0.024 (this would correspond 
to the "11 yr" sunspot cycle on the Sun). Then a massive 
alpha-perturbation was imposed in one longitudinal hemi- 
sphere, rotating with the angular velocity at the base of the 
convection zone (r = 0.7), for a time interval of length 0.005 
(i.e. about 0.2Pe). The perturbation was then switched off. 
The time evolution of the energies in the axisymmetric and 
nonaxisymmetric (overwhelmingly in mode m = 1) parts 
of the field are shown in Fig. [TJ 

Fig. [T] illustrates two relevant points. A nonaxisymmet- 
ric perturbation in alpha rapidly generates a (weak) nonax- 
isymmetric field. When the driver is removed (at time 0.605 
in Fig. Q]), the nonaxisymmetric field decays with a decay 
time of about 25% of the cycle period (i.e. of Pe). (Because 
of the disparate magnitudes of the quantities plotted in 
Fig. [T] the oscillation in the axisymmetric field is barely 
visible.) We do not pretend that this model is in any way 
realistic, but it does illustrate the points that a nonaxisym- 
mctic field is quickly established by the perturbation, and 
decays much more slowly. This is quite consistent with pre- 
vious studies which indicate that axisymmetric solutions 
for afl dynamos in spherical geometry are very stable to 
nonaxisymmetric perturbations, e.g. Radler (1986), Moss 
et al. (1991b). 

3. A scenario for reversal: mean magnetic field plus 
fluctuations 

A minimal extension of standard mean-field dynamo theory 
that includes solar magnetic reversals with a nonvanishing 
magnetic moment can be presented as follows. Mean-field 
dynamo theory assumes that, apart from the mean mag- 
netic field B which is considered to be the large-scale mag- 
netic field, magnetic field fluctuations b are also present, 
so the total magnetic field H = B + b. The mathematical 
expectation of b is zero but the spatial average of b re- 
mains finite. This is because the number of convective cells 
N, while large, is not so big that the spatial averaging of b 
(which scales as h/y/W) effectively vanishes . 

The number of convective cells N participating in dy- 
namo action in the convection zone can be estimated as fol- 
lows. Taking for estimates the supergranulation scale as 30 
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Fig. 1. The evolution of energy in the axisymmetric field 
(upper curve) and nonaxisymmetric field (lower curve). At 
dimensionless time 0.6, the steadily oscillating axisymmet- 
ric solution (period Pe ~ 0.024) is perturbed by the ad- 
dition of a nonaxisymmetric part to alpha between times 
0.6 and 0.605. The nonaxisymmetric perturbation is then 
removed, and the nonaxisymmetric field decays. 



Mm, we find that there are about f 4 supergranules at the 
solar surface or about (5 ~ 7) x f 4 in the whole convection 
zone. However supergranulation is a surface phenomenon 
whose spatial scale is probab ly controlled by the depth 
of the He ++ ionization region (Novemb er et al. 1981[) . The 
characteristic scales of the deeper convection are larger, so 
that TV sa 10 4 seems to be a reasonable estimate for the 
number of convection cells in the Sun. 

Knowing the magnetic field H we can estimate the 
dipole magnetic moment m as 



m 



const J Hd 3 x = m + amib / Vn . 



(1) 



where m is the magnetic moment of the mean field B, mi 
is a random vector of unit length, b is the r.m.s. value of 
the magnetic fluctuations and a is a numerical constant. 

Two mechanisms producing fluctuating magnetic fields 
are known: the small-scale dynamo and the wiggling 
of large-scale field lines by turbulent convection. The 
observation-based upper limit on the amplitude of the sur- 
face fields produced by the small-scale dynamo in the Sun 
is only about 3 G (Stenflo 2012 ). The contribution of such 
small fields to the global magnetic dipole moment of Eq. (JXJ) 
can be neglected in view of the large number N of con- 
vective cells. We therefore consider the distortion of global 
field lines by convection as the primary source of fluctuating 
fields. 

Assume that b is of order of the field strength of the 
toroidal magnetic field Bt- Such an estimate looks plausi- 
ble in the framework of the standard ideas of mean-field dy- 
namos and is supported by the analysis of unipolar sunspot 
group statistics (Khlystova & Sokoloff 2010; Sokoloff & 
Khlystova 2011). The poloidal (polar) fie ld B P of the Sun , 
which determines m, is ^1G (cf., e.g. |Hathaway 2010| ). 
The toroidal field near the base of convection zone is be- 
lieved to be at least 1000 times stronger. Near the sur- 
face, however, the toroidal field is probably weaker and 
the poloi dal field inside the convective zone is probably 
stronger (Kitchatinov fc Qlemskoy 2012). Taking for esti- 



mates Bp = 0.03-Bt and V N = 100 we find that the first 
term in the r.h.s. of Eq. (JXJ) to be several times larger than 
the second. In other words, it is the contribution of the 



mean magnetic field that determines the total solar mag- 
netic moment far from the instant of reversal. Of course, 
the direction of m does not coincide with that of m, i.e. the 
direction of the rotation axis precisely, but the angle 9 be- 
tween m and rotation axis is quite small, tan 6* = 0.1 ~ 0.2, 
i.e. 0«5°-lO°. 

The situation near the time of magnetic field reversal 
is quite different. Then m vanishes and the total magnetic 
moment m is determined by the second term in the r.h.s. 
of Eq. ([TJ). m becomes weak but does not vanish. With the 
above estimates, the magnetic moment at the instant of re- 
versal becomes several times weaker than its characteristic 
value far from the reversal epochs. 

When the Sun is still far from the instant of reversal, the 
first term in the r.h.s. of Eq. (1) is larger than the second. If 
mi is directed to the North pole, the direction of the total 
magnetic moment is aligned closely with this pole. For the 
same reason, m is aligned closely to the opposite pole after 
the reversal, when m is again large. The dipole's track from 
one pole to the other is determined by the direction mi, 
i.e. by the nonaxisymmetric part of the total magnetic field 
H. H is a physical field, and convection cannot destroy it 
immediately. This follows from the results of Sect. II where 
we find that the flow destroys a nonaxisymmetric field in 
a time of about 1-2 years, and is why the motion from one 
pole to the other is quasiregular. It is not fully regular just 
because the flow does not instantly destroy deviations from 
axial symmetry. 

In summary, m is determined by the sum of the dom- 
inant fluctuations (nonaxisymmetric) and the mean field 
contribution (axisymmetric). The nonaxisymmetric fluctu- 
ations individually decay on a timescale Si t r , and their 
sum is determined predominantly by the several most re- 
cent contributions. As the mean field component weakens, 
changes sign on passing through zero and then grows again, 
the sum m moves from being directed outwards in one hemi- 
sphere to pointing out from the other, on an irregular path 
essentially determined by the dominant fluctuation. 

The relative contribution of the second, fluctuating term 
in the r.h.s. of Eq. ([1]) can also be estimated as a ratio of 
reversal time to the cycle length. Taking 1-2 years for the 
reversal time t r and 11 years for the cycle length, we find 
again that the fluctuating part of m is up to an order of 
magnitude smaller than its mean part. 

According to Sect. II a nonaxisymmetric feature, and 
in particular mi, can survive for up to 1-2 years, i.e. the 
time t r required for a reversal. On the other hand, the con- 
vection turnover time t* estimated for the vortices of the 
largest scale is 1-2 solar rotation (cf. IStix 1991"]) . i.e. a few 
months - substantially less than the reversal time. Indeed, 
observers refer to the so-called 'fast' global changes on that 
time scale (see e.g. Hoeksema, 2006). Magnetic fluctuations, 
being determined by cumulative action of the velocity field, 
have a longer memory time, t r , than that of the convection 
itself (t*); however some rapid changes on the time scale t* 
have to be expected. 

The presence of two memory times, t r and t* , can ex- 
plain the irregular, chaotic track of the magnetic moment 
during the reversals. Observers (e.g. Livshits and Obridko 
2006) stress this feature of the reversals. 



3 



Moss et al.: Reversals of the solar dipole 



4. Discussion 

We have demonstrated that a minimum extension of the 
standard mean-field solar dynamo theory to allow for fluc- 
tuations is sufficient to describe magnetic field reversals 
with nonvanishing magnetic moment. With this in mind 
the main concepts of mean-field dynamos can be retained 
to explain the phenomenology of the solar cycle with rea- 
sonable accuracy; further development looks desirable. 

The initial formulation of mean-field dynamos included 
magnetic fluctuations at least implicitly. Hoyng (1987) 
stressed their role explicitly and Choudhuri (1992), Moss 
et al. (1992) and Hoyng (1993) applied this idea to expla- 
nations of solar cycle variability. The contemporary level of 
knowledge of fluctuations in the dynamo governing param- 
eter suggests that they may provide a scenario explaining 
solar Grand Minima (e.g. Moss et al. 2008, Usoskin et al. 
2009) and the Waldmaier relations for the time dependence 
of the solar cycle (Pipin et al. 2012). We conclude that the 
idea seems relevant to explaining magnetic reversals also. 

Note that a similar problem is present in the geody- 
namo. The geomagnetic dipole is currently inclined to the 
Earth's rotation axis at an angle 9 = 11° and, according 
to the paleomagnetic data (e.g. Laj et al., 1991) the geo- 
magnetic field has reversed from time to time in geological 
history, with nonvanishing magnetic dipole moment. It is 
reasonable to assume that, magnetic fluctuations underlie 
this phenomenon also. 

In summary, we feel that recognizing the importance of 
fluctuations in the context of the solar dynamo mechanism 
may provide a crucial step to explaining a range of non- 
periodic phenomena. 
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